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Abstract. This paper deals with a remarkable integrable discretization of the so(3) 
Euler top introduced by Hirota and Kimura. Such a discretization leads to an explicit 
map, whose integrability has been understood by finding two independent integrals of 
motion and a solution in terms of elliptic functions. Our goal is the construction of its 
Hamiltonian formulation. After giving a simplified and streamlined presentation of their 
results, we provide a bi-Hamiltonian structure for this discretization, thus proving its 
integrability in the standard Liouville- Arnold sense. 

1. Introduction 

This paper deals with a remarkable integrable discretization for one of the basic inte- 
grable systems, the three-dimensional Euler top, which describes the motion of the free 
rigid body with a fixed point. Equations of motion of the Euler top in the body frame 
read 

Xi = 0:1X2X3, X2 = 0:2X3X1, X3 = O3X1X2. (1) 

where x = (xi,X2,X3) G M'^, and the real coefficients Oj are parameters of the system. 
We will denote the vector of parameters by o = (01,02,03) G M^. Throughout this 
paper we will use an abbreviated notation, according to which (ijk) stands for any cyclic 
permutation of (123). Thus, system ([T]) takes with this notation the form 

Xi (y.iXjX]^. (2) 

The coordinates Xi stand either for the angular velocities in which case the coefficients 
Oj are given by 

ai = —J — , (3) 
or otherwise for the angular momenta Mj, in which case the coefficients Oj are given by 

ai = Y-T- (4) 

Here Jj are the principal moments of inertia of the body. The relation between the 
two formulations is given by Mj = /jfii. Integrability features of the Euler top include 
[1, 10, 11]: a bi-Hamiltonian structure, i.e. the existence of two compatible invariant 
Poisson structures on the phase space; two independent integrals of motion, which are 
in involution with respect to any of the invariant Poisson brackets; a Lax representation; 
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explicit solutions in terms of elliptic functions. For the reader's convenience, some of these 
features are briefly exposed in Sect. O 

The general problem of integrable discretization of integrable systems is dealt with 
in the monograph [11]. One finds there also a detailed exposition of an integrable dis- 
cretization of the Euler top, due to Veselov and Moser [8,12]. The basic feature of this 
discretization is that it comes from a discrete Lagrangian formulation on the Lie group 
5*0(3). Upon a reduction to so(3)*, it produces a correspondence, i.e. a multi- valued 
map, each branch of which is Poisson with respect to the Lie-Poisson bracket on so(3)*, 
like the original phase flow. Moreover, it shares the integrals of motion and the Lax rep- 
resentation with the original continuous time flow. This Lax representation is related to 
matrix factorizations. 

A class of discretizations of the Euler top sharing the integrals of motion with the 
continuous system has been introduced and studied in [2]. These discretizations are 
characterized by the equations of motion 



Here and below tilde denotes the shift t t + e in the discrete time eZ, where e is a 
(small) time step. In other words, in Eq. (and in similar situation throughout the 
paper) we consider Xi as functions on eZ, and we write Xi for Xi{ne) and Xi for Xi{ne + e), 

G Z. In Eq. ([5]), it is assumed that 7 ~ e/4 is some real- valued function on the 
phase space. Then the map ([5]) approximates, for small e, the time e shift along the 
trajectories of the continuous flow ([H). In [2] the functions 7 have been characterized for 
which the map {xi,X2,X3) 1—^ (xi,X2,X3) deflned by Eq. ([5]) shares the invariant Poisson 
structure with the continuous system. In particular, the function 7 for the Veselov- Moser 
discretization has been determined. A further integrable discretization of the Euler top 
belonging to the family was proposed in [4]. Interestingly, the simplest choice 7 = e/4 
leads to a a map which does not preserve the original Poisson structure. Discretizations 
([5]) share the Lax matrix with the continuous time Euler top. They are implicit, since 
these formulas represent a system of algebraic (nonlinear) equations for {xi,X2,xs) which 
does not possess a simple closed-form solution. 

The present paper deals with the following beautiful explicit discretization of equations 
of motion ([2]), introduced by Hirota and Kimura [5]: 



Here one can take 

5^ = f^, (7) 

we will adopt this choice for the vector of parameters S = (51,62,63) e throughout 
the paper. This discretization is explicit, since the algebraic equations ([6]) are linear 
with respect to {xi,X2,X3), and thus they can be solved in a closed form (see Sect. [3] 
for further details). Hirota and Kimura presented some of the integrability attributes for 
their discretization: two independent integrals of motion and a solution in terms of elliptic 
functions. Other attributes, like the Hamiltonian formulation and the Lax representation, 
has not been mentioned by them. The main goal of the present paper is to flll the 
first of these two gaps by providing a bi-Hamiltonian structure for the Hirota-Kimura 
discretization, and thereby to prove its integrability in the standard Liouville-Arnold 
sense. 



'^i -^i 'y^ii^'^j ~t~ '^j)('^k ~t~ -^k)' 



(5) 



■^i ■^i 6i(XjXif -\- XjX}^. 



(6) 
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We found it worthwhile to give also a simplified and streamlined presentation of the 
results found in [5]. Indeed, the discretization of the Lagrange top given by Kimura and 
Hirota later in [7], as well as some preliminary results by Ratiu [9], indicate that the map 
([6]) might be just a tip of an iceberg, a huge collection of discretizations of integrable 
systems of classical mechanics. We plan to develop this topic in a series of upcoming 
publications. 

It is an established fact that many of the most important integrable systems can be 
found in the classical literature on differential geometry. Usually this refers to solitonic 
partial differential equations, like the sine-Gordon equation, but it turns out to be true 
also for the integrable map ([6]): a 1951 paper in "Mathematische Nachrichten" by H. 
Jonas is devoted to a birational map {x,y,z) {x,y,z) given by 

X + X + yz + zy = 0, y + y + zx + xz = 0, z + z + xy + yx = 0, (8) 

which differs only unessentially from ([6]). The map (IHl) has an origin in the spherical 
geometry, (x, y, z) and (x, z) being the cosines of the side lengths of two spherical 
triangles with complementary angles. Jonas' results include integrals of the map (IHl) 
and its solution in terms of elliptic functions. Thus, [6] seems to be one of the earliest 
precursors of the theory of integrable maps. 

2. EULER TOP 

The aim of this Section is to recall some of the main features of the integrable continuous- 
time Hamiltonian flow 

Proposition 1. Let j3 = /32, /Ss) G &e a constant vector. A quadratic function 

HW = ^{f3,xl + /32xl + /3sxl) (9) 

is an integral of motion for ^ if and only if (3 ^- a, i.e. if /5iq;i + I320t2 + /^scts = 0. 
Proof: An easy computation based on Eq. ^ shows that 

^if(^) = + /32a2 + l32,a:i)xiX2X2,. 

□ 

Since the orthogonal complement of the vector a is two-dimensional, there are two 
independent integrals of motion. It is sometimes convenient to use a special basis of the 
orthogonal complement just mentioned, consisting of vectors with one vanishing compo- 
nent. 

Corollary 1. The three quadratic functions 

Gi = ]^{ajxl- akx]) (10) 

are integrals of motion for Of course, only two of them are (linearly) independent 
since aiGi + a2G2 + Oi^G^ = 0. 

Notice that any function H^^^ is a linear combination of the Gi's: 

aiH^f"^ = /3jGk - (3kG,. 
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In the angular velocities formulation, a basis of the orthogonal complement a-^ can be 
chosen consisting of (3^^^ = {Ii, 12,13) and = {I'l^I'^^lD- In the angular momenta 
formulation, a basis of consists of (3^^^ = {l/h, l/h, ^/h) and P^'^^ = (1, 1, 1). 

Proposition 2. Let /? ± a, and let 7 = (71,72,73) € M'^ satisfy 

Oli = Pjjk - Pklj, (11) 

so that 7 ± a. Then the system ^ is Hamiltonian with the Hamilton function H^^'> with 
respect to the Poisson bracket 

= 7fcXfc. (12) 

Proof: A direct verification: 

= iPjlk - (3klj)xjXk = atXjXk- 

□ 

Propositions [H and [2] show the bi-Hamiltonian property of the Euler top. Referring to 
the angular velocities the system has two Hamiltonian formulations: 

H = lihnl + hnj + i^nl) with {a, fij = -^fi^. 



and 



2 \ i ' ^ z ' /■/• 

H = -{I^fll + lifll + I^nl) with {Qi, Qj} = -^Qk- 
2 lilj 



Referring to the angular momenta, the system also has two Hamiltonian formulations: 

^ \ -'2 / 

and 

H = -{M^ + Ml + Ml) with {M„ = ^Mfc. 

3. HiROTA-KlMURA DISCRETIZATION OF THE EULER TOP 

We now turn to the study of the map Though the vector of parameters 5 is 

arbitrary, we will think of it as related to a as in Eq. ([7]). 

3.1. Integrals of motion. An explicit form of this map can be easily obtained. Con- 
sidering Eq. (E]) as a system of linear equations for the updated variables Xj, one finds 
immediately its solution: 

— 81X3 —5iX2 

1 -hxi 

-63X1 1 

Note also that, considering Eq. (I6l) as a system of linear equations for Xj, one finds the 
alternative formula 

1 (5iX3 5iX2^ 

82X3 1 52^1 

^83X2 83X1 1 
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We will use the notation 

/ 1 -5iX3 -61X2" 

A{x,5) = -62X3 1 -§2X1 
\-63X2 -S3X1 1 

so that the equations of the map can be written as 

X = A'-^i^x, 6)x = A{x, —S)x. 

Proposition 3. The quantities 

1 - SkSix'^j 
hSjxl 



1 - 5i5ixi 



are integrals of motion for the map Of course, there are only two independent integrals 
since F1F2F3 = 1. 

Proof: Equation Fi = Fi can be re-written as 

(1 - 5k5ix]){l - Si6jxl) = (1 - 6i6jxl){l - SkSiX^j), 

which is equivalent to 

6j{xl - xl) - 6k{x^j - x'^j) = 5i5j5k{x]xl - xpl), 

that is, to 

Sj{xk + Xk){xk - Xk) - Sk{xj + Xj){xj - Xj) = 6i6j6k{xjXk + XjXk){xkXj - XkXj). 
Using the equations of motion on both sides of the latter formula, we arrive at 

(^Xk ~t~ Xk^i^XiXj ~\- XiXj^ ~\~ Xj^(^XkXi ~\- XkXi^ (x^ Xi^(^XkXj XkXj^^ 

which is an algebraic identity. 

□ 

The relation between Fj's and the integrals of the continuous time Euler top is straight- 
forward: 

F, = l + ^G, + 0(e^). 

Corollary 2. Let /? ± 5. Then the following three functions are integrals of motion for 
the map 



1 - 6j6kxf' 



-'j^k-^i 

where the common numerator H^^"^ is an integral of the continuous time Euler top given 
m Eq. ^. 
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Proof: We show that iff ^ can be expressed in terms of the Fj's given in Eq. ffT^ : 

' ' ~ 1 - 6,6kxj 

PjjSjx'^j - 6jX^) + PkjSjxl - 6kx]) 
1 - 6j6kxf 

/3, l-h5.x] \ ^hL_l-h5,xl 



Sk \ 1 - 5jSkX^ J 5j \ 1 - 5j5kX^ 



□ 



3.2. Invariant volume form. Next, we estabhsh the existence of an invariant measure 
for the map 1^. Let us first give the following useful Lemma. 

Lemma 1. For the map ^ the following holds: 

Xi SiXjXk Xi ~t~ SiXjXk 



1 - 1 - 6j6kx'^ ' 

and, as a corollary, 

(^Xi SiXjXk^ {^Xi ~t~ S^XjXk^ 



(14) 



(15) 



(1 - 6i6kx]){l - 6i6jxl) (1 - Si5kx]){l - Si6jxl) ' 
Proof: We prove, for instance, Eq. ( fT4l) . It is equivalent to 

{xi - 6iXjXk){l - 6j6kxf) = {xi + 6iXjXk){l - 5j5kxl), 

or to 

Xi Xi SiXjXk SiXjXk S j5kXiXi{xi Xj) 5i5j5k{XiXjXk ~\~ x,^x jXk") ■ 

Upon using equations of motion ^ on both sides of the latter formula, we find that it is 
equivalent to 

{Xj Xj^iyXk Xk) ^ j^kiXiX j ~\- XiXj^iyXiXk ~\~ XiXk^^ 

which is a direct consequence of Eq. ([6]) . 

□ 

Now we are in the position to prove the following claim. 
Proposition 4. There holds: 

ox (p[x) 

where 0(x) is any of the functions 

<f){x) = {l-6i6jxl){l-6j6kx^,), (16) 

<P{x) = (1 - 6,6,xlf. (17) 

(The ratio of any two different functions (f){x) is an integral of motion for ^ due to 
Proposition^^ . Equivalently, the three-form 

Q = —rr dxi A dx2 A dx^ (18) 

0(X) 
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is invariant under the map 

Proof: First of all, we derive the following formula for the Jacobian of the map ([6]): 

dx _ detA{x,-6) 
"^^^dx ~ det A{x,S) • ^ 

Indeed, differentiating Eq. with respect to Xi,X2,X3, one obtains the columns of the 
matrix equation 

-61X3 -SiX2\ Q~ / 1 5lX3 5iX2^ 
-hXl 1 J ^ \S3X2 63X1 1 

Computing determinants leads to Eq. (|T^ . which can be written in length as 
, dx 1 - Sj6kX^ - 5i5kx] - Si6jxl + 26i6j6kXiXjXk 

Q-Gt = 

dx 1 - SjSkxf - SiSkx'j - 5i5jxl - 2SiSjSkXiXjXk 
_ (1 - 6i6kxJ){l - Si6jxl) - 6j6k{xi - diXjXkf 
(1 - 5i5kx'^){l - Si6jxl) - Sj6k{xi + diXjXkY' 

Now the claim of Proposition with as in Eq. ffTUl) , say, follows from Eq. ffT^ . 

□ 

3.3. Invariant Poisson structure. In the construction of an invariant Poisson structure 
for the map ^ we shall use the following results from [3] (Proposition 15 and Corollary 
16 there). 

Let / : M — i> M be a smooth mapping of an n-dimensional manifold M, and let VL be 
a volume form invariant under /, i.e., f*VL = Q. Define u to be the dual n- vector field to 
Q such that ujQ = 1. Here the symbol j denotes the contraction between multi-vector 
fields and forms. If h, ■ ■ ■ , In-2 are integrals of / with dli A ■ ■ ■ A dIn-2 7^ 0, then the 
bi- vector field a = uj dli ■ ■ -j dIn-2 is an invariant Poisson structure for /. If Ji, . . . , J„_2 
is another set of independent integrals and r = cuj dJi ■ ■ ■ j dJn-2 is the corresponding 
Poisson structure, then a and r are compatible, i.e., for any constants a, 6, the bi- vector 
field aa + 6r is a Poisson structure, again. 

In particular, for n = 3, if a three-form (|T8l) is invariant under /, so that the dual 
tri-vector field is given by 

, 9 d d 
uj = (l){x)— A — A — , 
OXi 0x2 0x3 

then for any integral / of / the bi-vector field 

,^ , / dl d d dl d d dl d d \ 

is an invariant Poisson structure for /, as well as any linear combination of such bi-vector 
fields. The Poisson brackets of coordinate functions are given by 

dl 

{xi,Xj} = ^{x)—. (21) 

Applying this result to the integrals logFi, logF2, logFs, with the three volume densities 
( TT6l) . we arrive at the following statement. 
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Proposition 5. The following brackets give compatible invariant polynomial Poisson 
structures for the map 



where Ci, C2, C3 are arbitrary constants. 

Notice that the Poisson brackets (!22l) yield three compatible polynomial Poisson struc- 
tures. Indeed, setting C2 = C3 = and Ci = 1, we get 

{xi,X2}i = 52X3(1 - ^s^i^s), {x2,a;3}i = 0, {x'i,xi}i = -5^X2{1 - 5i52xl), (23) 
setting Ci = C3 = and C2 = 1, we get 

{xi, X2}2 = -^1X3(1 - 525^x1), {x2, X3}2 = ^(l - 5i52xl), {^3, xja = 0, (24) 
while setting Ci = C2 = and C3 = 0, we get 

{a;i,a;2}3 = 0, {x2,X'i}'i = -52Xi{l - 5^5ixl), {x'i,xi}'i = 5iX2{l - 525^x1). (25) 

It is easy to verify that the brackets ( l23l) . ( l24l) . ( !25l) admit as Casimir functions the 
integrals Fi, F2, F3, respectively. 

In the continuous limit e — > these three brackets correspond to the invariant lin- 
ear brackets {■, ■j^'''^ of the Euler top, given in Proposition [21 with 7 = (0, —0:3, 0:2), 
7 = (^3,0,— tti), and 7 = (— a2,ai,0), respectively. Clearly, these three linear brack- 
ets are linearly dependent. On the contrary, the three polynomial brackets (!23l) . (|24l) . 
(1251) are linearly independent, if one considers linear combinations with scalar coefficients. 
However, they become linearly dependent, if one considers more general linear combina- 
tions. Indeed, the volume density (f) in Eq. (l20l) can be multiplied by an arbitrary integral 
without violating the Poisson property. Thus, in formulating the compatibility property 
of such Poisson tensors it is natural to consider their linear combinations with coefficients 
being integrals of motion rather than just numbers. In particular, the linear combination 
of the brackets (I2SD, (El, (ESD with the coefficients 



vanishes, so that there are only two independent brackets among them. 

3.4. Explicit solutions. Explicit solutions were given in [5], but it has not been ex- 
plained there how to determine the parameters of the elliptic functions involved in their 
formulas, using the initial conditions. We would like to fill in this gap here. We use the 
following addition formulas for the Jacobi elliptic functions: 



{xi, Xj} = Ci6jXk{l - 6k6iX^j) - Cj6iXk{l - Sj6kxf), 



(22) 




cn(^ + v) 



cn(^ 



v) 



2 sn ^ dn ^ sn 77 dn 77 



sn({ + v) 



sn(^ 



v) 



1 — k'^sn'^C, sn^r/ 
2 cn ^ dn ^ sn 77 



dn(^ + v) 



dn(e 



v) 



1 — /c^sn^^ sn^r/' 
2/c^ sn ^ cn ^ sn rj cn rj 



1 — fc^sn^^ 



ON THE HAMILTONIAN STRUCTURE OF HK DISCRETIZATION OF THE EULER TOP 9 

and the related formulas 

2 sn ,^ dn ^ cn r/ 



sn(,^ + ri)dn{^ ~ ^) + sn(^ — ?7)dn(,^ + rj) 
cn(^ + r?)dn(^ -v) + cn(^ - r?)dn(^ + r]) 
sn(^ + ri)cn{^ ~ v) + sn(^ — '^)cn(^ + rj) 



1 — k'^sn'^^ sn^T^' 

2 cn ^ dn ^ cn rj dn rj 
1 — k'^sn'^C, sn^r/ 

2 sn ,^ cn ^ dn r/ 
1 — k'^sn'^C, sn'^r] 



Assume that the coefficients Si are given by formulas with coming from Eqs. ([3]) or 
(Hj) with h < I2 < h, so that 

Si < 0, ^2 > 0, 53 < 0. 
Then the above addition formulas suggest to look for the solution in one of two forms: 
Xi = Aicn{iyn + ipo), X2 = sn(z/n + (/jq), X3 = A3 dn{iyn + (fo) , (26) 

or 

Xi = Aidii{i^n + (fo), X2 = A2sn{iyn + (po), X3 = ^3 cn(z/n + yjo), (27) 

with u being a parameter to be determined and ipo an arbitrary phase. Both possibilities 
fl^Uj) and fl?r|) are realized (in different regions of the phase space). Consider first the 
possibility fl^ . It is easy to see that equations of motion are satisfied by functions 
(!26|) . if and only if the following conditions hold [5]: 

,- . . cn(z//2) , , 

A,_ . M.^3 ""^^f,y^' . (29) 

sn(i//2) 

- -^-■^■■^- .^sn(t/2)on(./2) - 

The amplitudes Ai should be determined from the values of the integrals of motion. 
Substitute the ansatz fl2Ul) into the integrals ffT^ . then a direct computation based on the 
relations cn^,^ = 1 — sn^^ and dn^^ = 1 — A;^sn^^ leads to 

,42 _ 1 - ,2 _ 1 - -^3~' /.2 _ 1 - -^r' 

^2^)3 ^1(^3 

and 

,0 1 - 



Thus, this ansatz holds, if and only if Fi < F3 ^ < 1, that is, if F2 > 1. With the values 
just found, relations fEHjl - flHUl) lead to 



sn2(z//2) = 1 - Fi. 
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Turning to the possibility (1271) (omitted in [5]), we find that equations of motion 
are satisfied by functions ( 1271) . if and only if the following conditions hold: 

A K A A dn(t^/2) 

= fc^sn(./2)cn(./2) ' ^^'^ 

cn(z//2)dn(z//2) 

A2 = d2AiA3 — — , (32) 

sn(z//2) 

A X A A cn{u/2) 

A3 = -(33^1^2^——-—-. (33) 
sn(z//2)dn(z//2) 

Substituting the ansatz ( l27j) into the integrals (fT3l) . we find: 

42 _ 1 - 42 _ 1 - -^1 42 _ 1 - -^r' 

O2O3 0163 6162 

and 



1 - f3-' ■ 

Theferore, this ansatz holds, if and only if F-^^ < Fi < 1, that is, if F2 < 1, and then 
relations (!3Tl) -( |33l) lead to 

sn\u/2) = 1 - F^\ 

Thus, in both cases all parameters of the solution are expressed in terms of the initial 
data (more precisely, in terms of the integrals of motion). 

4. Concluding remarks 

In this paper, we studied a remarkable birational map of M^, which serves as an in- 
tegrable discretization of the Euler top, on one hand, and plays a role in the spherical 
geometry, on the other. Along with a streamlined presentation of results obtained previ- 
ously in [6] and in [5] , namely the conserved quantities and the solution in terms of elliptic 
functions, we found an invariant volume form and a family of compatible invariant Pois- 
son tensors for this map. Thus, it becomes a well-established representative of integrable 
maps, with a standard definition of integrability in the Liouville- Arnold sense. One more 
standard attribute of integrable systems remains to be found for this map, namely the 
Lax representation. This would provide a key to understanding the nature of analogous 
discretizations proposed in [7], [9], which seem to belong to the most mysterious objects 
in the universe of integrable maps. 
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